Rapidly-converging methods for the location of quantum critical points from finite-size data 
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We analyze in detail, beyond the usual scaling hypothesis, the finite-size convergence of static quantities 
toward the thermodynamic limit. In this way we are able to obtain sequences of pseudo-critical points which 
display a faster convergence rate as compared to currently used methods. The approaches are valid in any spatial 
dimension and for any value of the dynamic exponent. We demonstrate the effectiveness of our methods both 
analytically on the basis of the one dimensional XY model, and numerically considering c = 1 transitions 
occurring in non integrable spin models. In particular, we show that these general methods are able to locate 
precisely the onset of the Berezinskii-Kosterlitz-Thouless transition making only use of ground-state properties 
on relatively small systems. 

PACS numbers: 05.70.Fh, 64.60.an, 75.40.Mg 



i 

C 

O 

o 



> 
00 

m 
rn 

o 

oo 
O 



X 



I. INTRODUCTION 



In the study of physical properties of phase transitions a 
basic prerequisite is a reliable method to locate the critical 
point, whenever the latter is not known a priori from sym- 
metry or duality arguments. Typically, in numerical or even 
experimental studies on finite samples one obtains a sequence 
of pseudocritical points (in the sense specified below) to be 
extrapolated to the true critical point in the thermodynamic 
limit (TL). The extrapolation may be done with some polyno- 
mial fit in the inverse size of the system or, better, exploiting 
some fitting function derived on the basis of a scaling ansatz 
or through the renormalization group (RG). The point is es- 
pecially relevant in the context of quantum phase transitions 
(QPT) in lattice systems where the exponential growth of 
the dimension of the Hilbert space with the number of sites 
is a strong limitation on the accessible sizes with the current 
computational power and algorithms. One of the most used 
algorithms is still the Lanczos method for the virtually exact 
extraction of the low-lying energy levels; in the most favor- 
able case of spin- 1/2 models one cannot go beyond some tens 
of sites. This limit can be moved to maybe a few thousands of 
sites using the so-called density matrix renormalization group 
(DMRG) J2l that has become the method of choice for ID 
problems due to its high level of accuracy. Nonetheless, if one 
considers two or even three dimensional systems the situation 
is much worse: with the Lanczos algorithm the largest lat- 
tices have only a few sites of linear extension and the DMRG 
is not particularly efficient. At present, the only other choice 
is Quantum Monte Carlo (QMC) (see, for instance, JUl) that, 
however, suffers from a sign problem in the case of fermionic 
or frustrated systems and does not reach the level of accuracy 
of the DMRG. Very recently there have been attempts to ex- 
ploit both DMRG-like features and the QMC sampling tricks, 
to design hybrid methods [fl Ht] that are however still under 
verification. 

It is generally believed that a sequence of pseudocritical 
points, for example the loci of maxima of finite-size sus- 
ceptibilities, converges to the critical point as a power law 
L~ x , with a so-called shift exponent A given by the inverse 



of the correlation length exponent v. Hence, generally speak- 
ing, the larger is v the slower is the convergence. This 
difficulty reaches its maximum for Berezinskii-Kosterlitz- 
Thouless (BKT) transitions, in which the correlation length 
diverges with an essential singularity or, loosely speaking, 
"v = oo". However, already in the seminal paper by Fisher 
and Barber @|, it was pointed out that the relation A = \Jv is 
not always valid and A depends, among other factors, on the 
boundary conditions. 

The most used method to locate quantum critical points 
in d = 1 by means of finite-size data is the so-called Phe- 
nomenological Renormalization Group (PRG), reviewed for 
instance in 10]. Another convenient approach, the Finite-size 
Crossing Method (FSCM), was recently proposed in @]. The 
aim of this paper is to improve both of them by means of cri- 
teria that produce sequences of pseudocritical points that con- 
verge more rapidly. We will show that in our sequences the 
shift exponent will have the form A = Ao + Ai/V, where Ao 
and/or Ai are larger than the corresponding values in the usual 
methods and therefore allow for a better convergence. 

The paper is organized as follows. In Section [TT] we il- 
lustrate the general arguments leading to the enhanced se- 
quences, both in the framework of the FSCM and of the PRG 
(Subsec. Ill Dl >. Special cases as the BKT transition (Sub- 
sec. IIICI l and that of logarithmic divergences (Subsec. Ill Bl > 
are discussed separately. In Section [Til] we illustrate the use- 
fulness of the methods on the hand of analytic and numerical 
tests. In Subsec. UlI Al we treat the XY spin-1/2 chain, using 
a series of exact calculations reported in the Appendix. Then 
we move to two cases of spin chains for which no exact solu- 
tion is available: in Subsec. UlIBI we consider a spin-1 model 
with anisotropics in a parameter range that gives rise to a large 
value of v and in Subsec. UllCl we study the spin-1/2 model 
with next to nearest neighbor interactions that is known to un- 
dergo a BKT transition. In this case we find a value for the 
critical coupling in agreement with the accepted one, which 
was found using a model-specific investigation of the excited 
states [9]. Section HVlis devoted to conclusions. 
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II. DERIVATION OF RAPIDLY-CONVERGING 

SEQUENCES 

We consider systems in d spatial dimensions of linear size 
L and periodic boundary conditions (PBC). Let the transition 
be driven by a linear parameter g such that the Hamiltonian is 

H(g) = H + gV. 

Dealing with QPT we consider the case of strictly zero tem- 
perature, T = 0, even if the arguments presented below can 
be simply extended to the finite-temperature case, replacing 
the parameter g with T (and without using the dimensional 
crossover rule used below). The free energy density reduces 
to the ground-state (GS) energy density which, close to the 
critical point g c , shows a singularity in the second (or higher) 
derivatives with respect to g: 

£d(ft(s)> = e (#) = e rc g (3) + e sing (£(»), 

where £ m \g — g c \~ v = t~" is the correlation length. Note 
that, as a consequence of the scaling hypothesis, the singular 
part of the energy e S j ng is a universal quantity that depends 
only on £, the relevant length scale close to the critical point. 
Hence, e s j ng may be considered quite in general an even func- 
tion of (g — g c ) that vanishes at the critical point. 

On the other hand, the bulk energy density at the critical 
point behaves as (Privman-Fisher hypothesis) 

e(g c , L) = eoo (g c ) - L-^F{g c ), (1) 

where F(g) is a sort of Casimir-like term that may depend 
on the actual geometry of the lattice. Note that this hypoth- 
esis has to be changed properly if one or more of the spa- 
tial dimensions are of infinite extent. Moreover, Eq. (fTJ has 
been written in analogy with Eq. (11.29) of ref. IIIOII using 
the dimensional crossover rule according to which the parti- 
tion function and the thermodynamic (static) properties of a 
rf-dimensional quantum system are equivalent to those of a 
id + O-dimensional classical counterpart 01 OH, wner e C is 
the dynamic exponent QJ. Then, for the implementation of 
our methods we need to know by some other means the value 
of £ relating the energy gap A and the correlation length £: 
A oc £~^. Typically, but not always, energy and momentum 
in the continuum limit at the critical point satisfy a linear dis- 
persion relation, E = vk, for small k so that £ = 1 and a 
relativistic effective field theory can be used to describe the 
universal features of the transition. In d = 1 the scale invari- 
ance at the critical point is often sufficient to imply also con- 
formal invariance (see Ch. 2 of frill ), thanks to which several 
exact results can be obtained using the powerful predictions 
of conformal field theories (CFT). For example, by mapping 
the space-time complex plane onto a cylinder whose circum- 
ference represents the finite chain of length L we can identify 
F(g c ) = ircv(g c )/6 in Eq. (Q]i where c is the central charge 
of the theory. In the RG sense, moving away from critical- 
ity corresponds to perturbing the CFT with a relevant operator 



that destroys conformal invariance. However, this is not the 
only effect of varying the microscopic parameter g out of g c : 
in general also the speed of elementary excitations gets renor- 
malized in the unperturbed CFT part. For this reason we say 
that v (henceforth F) depends on g, in the vicinity of g c . 

Scaling and dimensional arguments imply that, in the ther- 
modynamic, off-critical regime L 3> £, the singular part of 
the energy behaves as 

^sing ~ t (2) 

with a = 2 — [d + £) v. For a second order phase transition 
a < 1. 

After introducing the scaling variable z = {L/S,) 1 ^ = 
tL 1 ^, the finite-size scaling (FSS) theory asserts 0] that in 
a system of length L, 

e sine ^ CoL-^^ e (z) + ■ ■ ■ , (3) 

where Q e {z) is a universal function that, in the off-critical 
regime z 1, must behave as 3> e (z) ~ Z ( d +Qv m order 
to recover Eq. (f2]i . Instead, for L <C £ we are in the critical 
regime and $ e (z) behaves as an analytic function that van- 
ishes for z — ► 0. Here we assume that the leading term in 
<& e (z) is quadratic in z (see below why it cannot be linear), 
but the following arguments are easily generalizable to higher 
integer powers. The constant term Co$ e (0)/L +< * is already 
adsorbed in the nonuniversal part of the energy density at g c 
as shown in Eq. (Q]). 

Differentiating e(g) with respect to g, gives the mean value 
b = (V) / L d , whose singular part 6 s ; ng behaves as 

6 sing « sgn(.g - gjt 1 -* = sgn(.g - g^t^"- 1 . (4) 

Considering FSS for the combination of Eqs. (Q]i and (O and 
then differentiating we find 

b(g, L) = 600,^(3) + sgn(.g - g c )C L 1 '''-< d +O& e (z) 
_ L -(d+0 [F ' {gc ) + F"{g c ){g -g c ) + .~ }+0 

1 

where the subscript "00, reg" hereafter means regular in the 
TL. In order to write down the expression above we used 
d g = L x l v d z and assumed that the powers neglected in the 
last term are just larger than (d + Q. To illustrate this point 
we could consider the irrelevant operator with the smallest 
scaling dimension di rl . At first order in perturbation theory 
with the renormalized coupling g- m (L) = g- m {Q)L d+< ^~ di " 
the corrections to the GS energy density are of the form 
Ci„ghr(0)L~ di " so that e = d hl — d — £ > 0. Note also 
that the amplitude C m can vanish and 0(gf IT ) terms have to 
be included. For these and more details we leave the reader to 
Ref. lfl2ll . Generically we admit corrections with e > that 
may come either from irrelevant operator in the continuum 
theory or from lattice effects. The case e = corresponds to 
marginal perturbations and typically leads to logarithmic cor- 
rections. Notice that now the leading term of Q' e (z) is linear 
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in z in the critical region. If we would have admitted a linear 
term in $> e (z) then ^(O) ^ and a finite jump discontinuity 
at finite L would be present in b(g, L). 

We can also calculate a = (Ho}/L d = e(g)—g[de(g)/dg], 
yielding to a singular part that is similar to Eq. but with a 
changed sign. In fact, for g c ^ the leading singular parts 
of a and b must cancel in the sum that gives back the energy 
e(g) = a(g) + gb(g), which does not contain that singularity. 
In particular, the scaling is 

a(g, L) = aoo>Tes {g) - sgn(.g - g^gC^-^^z) 

-L-( d+ V [F{g) _ gF ' { g c) + (g - g B )] + ... 

The FSCM H identifies the critical point with the limit of the 
sequence g* L of single crossing points 

b{g* L ,L) = b(gt,L') 

with L' = L + SL. Applying this criterion to Eq. ©, we 
obtain (for SL «C L) 



9l - 9c 



(d + C)F'(g c ) _ L _ 2/ ^ ^ 2 



2C [§-(d + 0] 



This equation defines the shift exponent Afscm = and 
may converge very slowly when v 1. The extremely diffi- 
cult case is the BKT transition where formally v = oo, but this 
latter situation must be treated in a different way (see Subsec. 
HTCl 

Now, consider the quantity T(g, L, 7) = je(g, L) — b(g, L) 
and suppose to be able to tune 7 exactly at 



7 



F'(g c ) 
F(g c ) 



(6) 



It easily seen that T(g c , L, 7*) does not contain the Casimir- 
like term responsible for the critical point shift. In fact, the 
scaling of T is 

r(g,i )7 ) = r oo , rog (. g ,7)-sgn( 5 - gc )C L 1 ^- (d+c) ^(z) 
-L-^ [ 7 F(. 9 ) - F'(g c ) + 0(g - g c )}+0 (V< W e >) . 

When 7 is equal to 7* given in Eq. ©, the critical point found 
by this crossing method is approached as g* L — g c ks L~ Afaat , 
with a shift exponent Af as t = 2/v + e = ApscM + e (again this 
holds true provided that v ^ 2/(d + ()). The additional term 
e allows, in general, for a better convergence of the sequence 

A possible algorithm for finding numerically the critical 
point in such a way is the following. If g* L is at a crossing 
point, of r(<7, L, 7), then we have 



7(Sl,£) = 



b(gi,L + 6L)-b(gj,L) 
e{gl,L + 5L)-e(g* L ,Ly 



Notice that putting 7 = is equivalent to the FSCM applied 
to b(g, L) and the denominator has definite sign about the crit- 
ical point. Now we find g* L requiring that 7(^2, L — 5L) = 



b(gj,L) - b(gj,L - SL) = b(gj,L + SL) - b(gj,L) 
e(g* L , L) - e{g* L , L - 5L) e(g* L ,L + SL) - e(g* L , L) ' 



or in the continuum version 

d L b(g,L) 



_ d 2 L b{9,L) 
d L e(g,L) d 2 L e{g,L) 



(7) 



(8) 



Calling t = (g — g c ) = tsgn(g — g c ) and rewriting only the 
essential terms in the scaling ansatze we have the simplified 
forms 

e(g, L) = eoo . Jce (g) - L-( d +OF(g) + D 1 L~^+^, 
b(g, L) = b^ rcg (g) + L-< d +« [2C i 2/,y r - F'(g) 

+D 2 L~ (d+ t +t \ 
Now, putting these two relations in ^ we obtain 



9l-9c 



D 2 F(g c )-D 1 F'(g c ) ve(d + ( + e) w 
4C F(g c ) 



(9) 



that gives a shift exponent Af as t as anticipated. 

The main result of this section is the crossing criterion (0 
that identifies the rapidly converging sequence (0 to the criti- 
cal value. 



A. Homogeneity condition 

In the previous section we have shown how to obtain a se- 
quence of pseudocritical points with an improved shift expo- 
nent Afast ■ Here we provide another, yet simpler equation for 
the determination of the critical point. The resulting pseudo- 
critical sequence is characterized by the same shift exponent 
Afast ■ However in this case we are able to prove convergence 
toward the critical point even in the extreme case of a BKT 
transition (see Subsec. ITTOb . 

The idea is to require that at the critical point the L- 
dependent part of b is dominated by the Casimir-like term with 
power (d + () (see Eq. (0). This condition is translated into 
the requirement that the i-derivative of b is a homogeneous 
function of degree (d + £ + 1), i.e. 

(d + C + l)d L b(g* L , L) + Ld 2 L b(g* L ,L) = 0. (10) 

Consequently, the corresponding sequence of pseudocritical 
points g* L scales as 

_ eD 2 (d + C + e) . t 



9l 



(ID 



2Co|(d + C-|) 

The equation ( TTOb represents the homogeneity condition 
method (HCM) that we are proposing for the efficient loca- 
tion of critical points. We stress here that, being b (g, L) a GS 
property, this criterion does not require knowledge of excited 
states as in the case of the PRG method. This is a point in fa- 
vor to the HCM since excited states are typically assessed with 
less numerical accuracy. In addition, the HCM is superior to 
the PRG in that it produces a faster converging sequence (see 
Subsec. HTdTi. 
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B. Case v = 2/(d + Q with logarithmic divergences 

For completeness we consider the case v = that was 
excluded in the previous treatment. In this situation, the 
ansatz requires the inclusion of logarithmic corrections 

b(g,L) = 6 00 , reg (5)+C 1 \nLT-L-^F'(g)+D 2 L- ( - d +^\ 

The calculation of the critical point with the FSCM gives in 
this case 



whereas the HCM Eq. ([Toj yields 

eD 2 (d+C + e) L _ (d+(+e) 



9l - 9c 



Ci(d + Q 



These results are compatible with the exact calculations of the 
XY model (see Subsec. IIII Ab . Note that in general, we should 
to perform a similar calculation for 



P 



d + C 



peN + l, 



namely when the p-th derivative of the free energy diverges 
logarithmically. 



C. A scaling ansatz for the BKT case 

For d = 1 at the BKT the correlation function in the TL be- 
haves like £ ps cxp(a£ _cr ). In the typical example of the clas- 
sical two dimensional XY model it is known that a = 1/2. 
Instead, for the quantum Heisenberg model with frustration 
(which we will consider in Subsec. IIII Cb Haldane suggested 
a = 1 [9]. We also set £ = 1 because the (effective) di- 
mensionality in the BKT scenario is two. The singular part of 
the finite-size energy density now is conveniently expressed 
in terms of y = L/£ fjv|] so that 



; (g)+L- 2 [C * e (y) - F{g)]+0 (l^+< 



where $ e (?/) is a universal function that, in the off-critical 
regime y ^> 1, must behave as $ e (y) ~ y 1 ■ Again, in the 
quasicritical regime y <C 1 at any finite L the energy density 
and its derivatives must be analytic in (g — g c ). The value 
$ e (0) can be absorbed in F(g) and it can be checked directly 
that the first contribution has to be at least quadratic in t = 
\g — g c \ because otherwise a finite-size discontinuity in b(g, L) 
would be generated. For y <C 1 we adopt the following ansatz 
(justified from perturbed conformal field theory lfl2"tl ): 



e{g,L) 



--oc,rcg V 



(g)+L- 2 \K(at-° - \nLy n '° - F(g)] +. 



with K a constant and n an integer larger than 1 (n = 3 from 
Eq (22) in |H). Hence 



^oo,reg(<?) 




nKa 








r 


a 





0(L-< 2+6 >) (12) 



F'(g) 



Now we want to get rid of all the 0(L~ 2 ) contributions 
that "hinder" the location of finite-size pseudocritical points. 
Hence we first differentiate with respect to L to eliminate 
boc,rcg{g), then multiply by L 3 to isolate the term in square 
brackets in Eq.(fT2b and finally set to zero a further difference 
in L in order to drop F'(g). Formally, in the region L/t; <C 1, 
we can write down the condition 



d L [L 3 d L b(g* L ,L)] =0. 



(13) 



It is worth noticing that the latter condition is equal to the 
HCM Eq. ( fTOb when d = ( = 1. Treating L as a continuous 
variable one can read off the shift exponent for the sequence 
g* L that turns out to be A^ t T = e/(rc-l + cr);ife = then the 
corrections to scaling are also governed by (another) marginal 
operator and we expect g* L —g c ~ (mL)~( m+1 )/("~ 1+fT ) with 
m a positive integer (m = 4 from Eq. (22) in lfl2ll ). 

Dealing with numerical simulations it is very important to 
specify how one implements the finite-size differences in L. 
In fact, there are several finite-difference expressions used in 
the literature to express the derivatives and here the require- 
ment is that they all reproduce Eq. dT3l > in the limit L — > oo. 
For example, if one takes a uniform step SL then the following 
symmetric expression can be built 



b'(g,L) = 
b"(g,L) 



b(g,L + 6L)-b( gi L~6L) 
2SL 

b(g, L + SL)- 2b(g, L) + b(g, L - 5L) 



and then the precise condition to cancel the term 0(L 2 ) be- 
comes 

L 3 b"(gl,L) + [3L 2 - [5Lf] b'{g* Ll L) = 0. (14) 

In the limit of large L we recover ( fT3l as required. Clearly, 
the correct discretization prescription must be identified not 
only for the BKT; when d and/or ( are not one, the suitable 
variant of ( TT4b for finite L has to be adopted, with b" weighted 
by L^+C+i) and b' weighted by a polynomial in L of degree 
(d + () whose coefficients depend on (5L). 



D. PRG revisited 

The PRG method identifies the critical point with the limit 
of the sequence g* L of crossing points satisfying 



L<A (g* L ,L) - (L + 8L^A(g* L , L + 5L) = 0. 



(15) 
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Here A(g,L) is the finite-size energy gap of the spectrum for 
which we may adopt the following form 

A(<7, L) = L d [ ei (g, L) - e(g, L)] = L~<*{z) + L~^C X 

(16) 

with $(z) ~ $ + z$i + z 2 &2 + • • ■ and C x a prefactor 
depending on the excited state \x) we are considering. The 
standard PRG approach in d = 1 relies upon Eq. (15[ with 
C = 1. More generally, a first test to identify ( is done by 
plotting the usual scaled gaps LA and see if at the critical 
point they settle to a constant or not. If they do not, one is 
led to search for a better value of ( ^ 1 and solve (Q3) with 
the correct value of the dynamic exponent. Note that the cor- 
rections come from the irrelevant/lattice contributions and the 
Casimir-like term for the gap is actually the constant term in 
$(z). For quantum systems in d = 1, CFT (£ = 1) pre- 
dicts $o = 2 7T d x v(g), where d x is the scaling dimension of 
the operator \ that generates the excited state \\). In addition 
$i, $2 7 ■ ■ ■ can be computed in the framework of perturbed 
CFT, in the regime z <C 1 when a relevant operator (g — g^fR, 
with scaling dimension x = 2 — 1/v is added to the critical 
field theory. For instance, from Eqs. (7) and (10) of Ref. lfl2tl 
we have $i = b x /(2ir) x where b x is the structure constant 
that appears as prefactor of the three points correlation func- 
tion (ft(ri)x(r 2 )x(r 3 )). 

Now, in general the PRG method gives \g* L — g c \ f=a £~ a prg 
with Aprg = 1/v + e 111 ill - However, it may happen that 
$i = in the scaling function (examples are discussed in 
Subsec. IIIIBI and at the end of Subsec. IIII A I for a specific 
excited state of the XY model). The shift exponent in such a 
case decreases to Aprq — 1/v + e/2. Nonetheless there is a 
way to improve this behavior. In fact, the "extremum" (instead 
of the zeroes) of the quantity in left side of Eq.dTBTl with the 
ansatz dT&b is located exactly at the critical point z = 0. In 
order to appreciate a shift from criticality we have to include 
higher orders of r = (g — g c ) in the irrelevant non-scaling 
term, i.e. 

L C A( 5 , L) = $ + £ 2/ 'V$2 + (C x (0) + C' x (0)t) L-\ 

In this case, the convergence is \g* L — g c \ « L _Acx " with 
A ox t r = 2/v + e, which is better than usual PRG not only 
because of the double exponent but also thanks to the coef- 
ficient in front: it is proportional to C' x (0) which is usually 
a small quantity. It is also worth noticing that in this case 

Acxtr = Af as t. 



III. TESTING THE METHODS 

A. XY model in transverse field 

To check our methods analytically, we consider the ID 
spin- 1/2 XY model given by 



H = - 



L r 

£ 

3=1 



(1 + ??) 



V V 



Throughout the paper we will consider L even and PBC. This 
model can be solved exactly ifTH [l3"tl by means of a Jordan- 
Wigner transformation from spins to spinless fermions fol- 
lowed by a Bogoliubov transformation to arrive at a Hamil- 
tonian of free quasiparticles. The number of original fermions 
jV = J^i c l c i> i s n °t a conserved quantity, but its parity 
P = cxp(iirj\f) corresponds to a it -rotation around the z-axis, 
and therefore is conserved. One should beware of a delicate 
issue concerning boundary conditions. Starting with PBC in 
Eq. ( fTTb . the fermionic Hamiltonian turns to have PBC in the 
sector of odd parity P = — 1. Instead in the even parity sec- 
tor P = 1 - which comprises the ground state - antiperiodic 
boundary conditions must be used. In this sector the model 
becomes 



(18) 



where (3^ are Bogoliubov quasiparticles, k ranges in the first 
Brillouin zone and the dispersion relation is 



A(jfc) = 2yV + h 2 + (1 -?7 2 )cos 2 /fc + 2/icosfc. (19) 

For r\ 7^ this model displays an Ising transition at h = 1 
with exponents v = £ = 1, which means that the scaling 
variable is z = L(h—1). To test our ideas we need to 
calculate, for finite L, the GS energy (h, rj) and the aver- 
age potential (h,ij) = — {cr|} = de^/dh, in the quasi- 
critical region given by z <C 1. According to Eq. ( TT~8b 
the GS is given by e L (h,r]) = -l/2Lj2 kn A witri 
k n = (2n + l)ir/L, q = 0, . . . , L — 1. We then expand 
the argument of the sum up to the desired order in z. For our 
purposes we need eL(h,r]) up to 0(z 2 ) and bh (h,r/) up to 
O (z). The resulting sums are then evaluated with the aid of 
Euler-Maclaurin formula flijl . The results and the details are 
given in the Appendix[A] From Eqs. ( IA2t and JA31) . one sees 
explicitly that terms of order L~ 3 are absent both from e_L and 
b^ so that according to our definition Eq. (O we find e = 2. In 
passing we notice that, from Eq. iA2i . the Casimir-like term 
is — L -2 7r|?7|/6 consistent with the CFT formula —L~ 2 ttcv/6. 
In fact the central charge is c = 1/2 and, from the dispersion 
relation (JT9J the spin velocity turns out to be v = 2\r/\. 

Now, using Eqs. ( IA2t and ( IA3I ). we have all the elements 
to derive the sequences of pseudocritical points analytically. 
As far as the FSCM is concerned, the pseudo critical points 
are obtained imposing bi (h* L ) = b^+i {h* L ), or formally 
dibi (h%) = 0. Up to leading order in L the solution is 



h* L = l 



TT 

MJ 2 



(17) 



as already obtained in Ref. |§] for the Ising model (|^| = 1). 
This shows explicitly that Afscm = 2 consistent with our 
prediction Afscm — 2/v (remind that v = 1). 

For what concerns the "balancing trick" discussed in 
Sec.|II] we can show that the solution of Eq. (O, is given by 
7* = — 1/2?/ 2 . The pseudocritical points are then given im- 
posing diX{h* L ,L, 7*) = 0. In this case, at leading order, the 
solution is 

7tt 4 (27 ? 4 + , ? 2 - 3) r _ 4 



1 



720?y 2 



6 



which means Af ast = 4 once again consistent with our predic- 
tion Af as t = 2/v + e (remind that e = 2). From a numeri- 
cal point of view it is more profitable to use the HCM. From 
3d L b + Ld\b = we obtain 



B. c = 1 transition (non BKT) 

We choose the spin-1 J z — D model on a chain (d = 1) 
with PBC 



K. = 1 + 



7tt 4 (2^-3)^ 



720?y 2 



with the same exponent Af ast = 4. 

Now we proceed to discuss the PRG method for which 
knowledge of the lowest gap is required. The first excited 
state belongs to the sector with odd parity. Correspondingly, 
the finite size gap is given, besides a constant term, by the 
difference between two Riemann sums where the sampling is 
taken over odd and even (in units of ir/L) wavenumbers re- 
spectively. In analogy with Ref. lfl5ll the lowest gap can be 
eventually written as A L = 2(h - 1) + [T(L) - 2T(L/2)], 

where T(L) = (1/2) Y^Lq 1 A(ji/i). Again we refer to the 
Appendix for the details. Using the form of the gap Eq. (IA4b 
we can evaluate the various terms of the scaling function in 
Eq. (O 



7T|7?| 
2 



TT 3 [3 (h - 1) - 2 (2 + h) if + 8r/ 4 ] 



192 



V 



Hence <f>o 



7T|7J| 



2irvd x , consistently with Ref. 



for |?/| = 1 and with the known result d x = 1/8 for the 



first excited state of the c = 1/2 minimal model 111 111 . More- 
over, since $i ^ 0, according to the general analysis in Sub- 
sec. Ill Dl there is no advantage in searching the minima or 
maxima in r of the left hand side of Eq. ( fT3l l. Looking for 
the zeros by imposing the PRG equation Ol (LAl (h^)) = 0, 
yields 



( V - 3) 
48 M 



LT 



0(L~ 



This result for generic 77 7^ 0, shows explicitly that the PRG 
shift exponent is Aprg = 3 and proves a conjecture put for- 
ward by Hamer and Barber in Ref. (TBjj for the Ising model. 
From Eq. (|20| | together with Eq. ( TToT l we read e = 2 as pre- 
viously, so that the calculated exponent is consistent with the 
prediction Aprg = l/v + e. 

The explicit calculation of the various shift exponents, con- 
firms that the HCM given by Eq. (TTOb is superior to the 
standard PRG method. However, the things are different if 
we considered the excitation obtained creating a well-defined 
quasiparticle with the smallest momentum k = 2tt/L, instead 
of the very first excitation gap. In that case, the energy gap 
is given by Eq. $1% . Using the method of the extrema of the 
PRG quantity in the left hand side of Eq. ([L5l l, we would have 
obtained a convergence to the critical point as L~ 4 , while with 
the standard PRG only as L~ 2 . 



H = J^2§ r S j+1 + (J z - J)S*S* +1 + D(S z 3 f (21) 
3=1 

because the transition from the Haldane phase to the phase 
at large D is described by a c = 1 CFT (£ = 1) with con- 
tinuously varying exponents. The Hamiltonian (f2Tb has been 
used to describe the magnetic properties of different quasi ID 
compounds (see lfl6ll for a brief account). Here, by fixing 
J = 1 and J z = 0.5, for which it has been already esti- 
mated v = 2.38 lfr7ll . we wish to test the methods described 
above in a severe case in which a 1/y scaling would give sub- 
linear convergence in L. Indeed for the FSCM we expect 
Afscm = 0.84. The previous estimate using the log-log plots 
of the finite-size gaps was D c = 0.65 ifnll . while in Ref. lfl8tl 
it is found D c = 0.635 using the method of level crossing 
with antiperiodic boundary conditions, which is however spe- 
cific to this transition. The first irrelevant operator allowed 
by the lattice symmetries has scaling dimension K + 2 where 
K = 2 - i so that e = mm(K, eiattice) where with entice 
we denote the smallest exponent of corrections arising from 
lattice contributions at finite L that are not captured in the 
framework of the (relativistic) continuum theory. 

As discussed in Refs. lfl6l [l9tl this is a case in which the 
linear term $1 in the scaling function of the PRG vanishes 
because b x = for the sine-Gordon model, the effective field 
theory that describe the surroundings of the c = 1 line. So 
it is convenient to use also our improved version of the PRG 
as discussed in Subsec. Ill Dl The expected shift exponent is 
Aoxtr = 2/i/ + e. 

In order to have an idea of the range of values of L to be 
used let us imagine that e = K so that Af as t = A ox t r = 
2 + 1/u = 2.42. With L > 10 this exponent leads to vari- 
ations in the pseudocritical points D* L smaller than O(10 -3 ). 
Hence we prefer to illustrate the method with virtually ex- 
act numerical data obtained with the Lanczos algorithm using 
L = 8, 10, 12, 14, 16. The reason is that the DMRG would 
give rather accurate values for the energies but the estimates 
for b = ((S z ) 2 ) could not be sufficiently precise to appreci- 
ate the variations in Dl obtained from the crossings. So we 
use the DMRG only to extend the data to L = 18, 20 with 3 7 
optimized states. In any case the GS belongs to the S^ ot = 
sector. 

The FSCM with b = ((Sp 2 ) yields D^ SCM = 0.647 with 
Afscm = 0.79 in reasonable agreement with the CFT expec- 
tations. The high-precision procedure based on Eq. ^} yields 
D f c aBt = 0.633 with A fast = 7.6. The value of the shift ex- 
ponent is definitely larger than what expected, which could be 
due to the vanishing of the coefficient of the first irrelevant 
contribution with scaling dimensions K + 2 < entice- In any 
case, from Fig. Q]one can clearly appreciate that the sequence 
£)FSCM conver g es m0 re slowly than D f £ Bt . As far as the PRG 
is concerned, with the standard procedure of finding the ze- 
roes of Eq. O we get D™ G = 0.640 with A PRG = 1.50. 
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Figure 1 : (color online) The sequences of pseudocritical points ob- 
tained with the FSCM, the fast-convergent method of Eq. ((7} and 
the HCM described by Eq. d 14b . Moreover the sequences have been 
calculated with PRG methods, both standard and in our improved 
version. Continuous lines are algebraic best-fits to the data. 



Instead, using the improved method estimation by looking for 
the extremal value of the quantity on the left side of ( fl~5T > we 
find a rapidly converging sequence that however also oscil- 
lates between 0.636 and 0.635. Finally, a small oscillation 
in the fourth decimal place about 0.6305 is seen also in the 
sequence D^ om obtained through the homogeneity criterion 
O. As expected the latter sequence converges in a fashion 
similar to Z?^ st . All these results are summarized in Fig. [TJ 

The data show that subtracting out the terms that induce 
a slow convergence of the pseudocritical points, one remains 
with quantities that inherit a residual (possibly oscillating) in- 
dependence from the specific lattice model and that could be 
very difficult to account for. Other factors that affect the ex- 
trapolation of the critical points at this level of accuracy are 
the sampling 5D and the trade-off between computational ac- 
curacy and the maximum available size. 

To summarize we observe that our improved methods yield 
very fast convergence to the critical point that we estimate to 
be D c = 0.633 ± 0.02, consistently with Ref. CJ]. 



C. c = 1 BKT transition 

We consider now the spin- 1/2 Heisenberg model with frus- 
tration due to next-to-nearest neighbors interaction 

L 

H = ^ ■ Sj+i + J 2 Sj ■ Sj +2 . 

The model is equivalent to a 2-legs zigzag ladder with L/2 
rungs. The best estimate of the critical point was made 
by Okamoto and Nomura using a model-specific crossing 
method 0]; with exact diagonalizations up to L = 24 they 
determined J 2c /Ji = 0.2411 ± 0.0001. The model is gap- 
less for J2 < J2C and has a doubly degenerate GS in the TL 



for J 2 > J2c- Again the GS has S^ ot — 0. Without exploit- 
ing a priori information about the BKT character of transition 
(if not the value of the dynamic exponent £ = 1), we tested 
the homogeneity criterion (fT4l using ladders with PBC of up 
to L/2 = 16 rungs, an effective SL = 4 and 1024 DMRG 
states that ensure an accuracy of O(10~ 7 ) on the values of 
b = {Sj ■ Sj+2)- The results reported in Fig. [2] are encour- 
aging: while with the FSCM we would get no crossings at 
all, the zeroes of Eq. ( fT4l i yield a sequence of points con- 
verging to J 2 ~ 0.25. The main problem comes from the 
left side of the transition where the truncation DMRG error 
for L — 28 and 32 induces some oscillations on the plotted 
quantity. We content ourselves with linear fits in 1/L, If we 
exclude the point with L = 12 the fit is better even if we 
find J 2c /Ji = 0.2553 ± 0.0008; by selecting all the avail- 
able points, instead, the fit is visibly worse but the extrapo- 
lated value is Jo c / 1 J\ = 0.242 ± 0.006. in agreement with 
Ref. H]. As above, apart from the details of the extrapolation 
procedure, we see that the homogeneity criterion provides a 
viable procedure to locate the critical point in a BKT transi- 
tion, where almost all existing generic methods fail. We re- 
mark that this analysis is based solely on b(g, L), namely an 
observable evaluated on the GS, without invoking further as- 
sumptions on the nature of the excitations. 



IV. CONCLUSIONS 

Making only use of finite-size quantities related to the 
ground state, we show how to generate sequences of pseu- 
docritical points that converge very fast to the infinite-size 
critical point. The convergence is of the form L~ x with a 
shift exponent A. In this article we propose a homogeneity 
condition method (HCM) which is faster than the standard 
phenomenological renormalization group (PRG) in locating 
the critical points. Moreover its validity is more general as it 
can be applied without modification to the difficult case of a 
Berezinskii-Kosterlitz-Thouless transition. The homogeneity 
method requires only the knowledge of b(g,L) = (V)/L d , 
that is, the expectation value of the term that drives the transi- 
tion. 

We also presented an improvement to the PRG method, 
which allows, under certain conditions, to obtain pseudocriti- 
cal sequences characterized by the same shift exponent as for 
the HCM. However, this modification, relying on a particular 
form of the gap scaling function, is not valid in general. It 
holds true, for instance for the sine-Gordon model that under- 
lies a variety of transitions in (1 + 1) dimension. 

The formulations of the approaches are sufficiently general 
to be applied in any spatial dimensionality. Even if we are 
primarily interested in quantum phase transitions ideally at 
zero temperature, in principle the methods can be extended 
to problems of finite-temperature statistical mechanics. At 
variance with other accelerating methods found in the liter- 
ature, e.g. the van den Broeck-Schwartz or the Bulirsch-Stoer 
ones (reviewed in Ch. 9 of ifTlll ). the procedures presented 
here rely on the scaling behavior of thermodynamical quanti- 
ties expected from physical and renormalization arguments. 
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Figure 2: (color online ) Upper panel: searching for the zeroes of 
L 3 b" + (3L 2 - 16)6' according to Eq. (H with Ji = 1 and 5L = 4. 
Lower panel: extrapolations of the pseudocritical points J2,l (cyan 
line for all the points, blue line with L = 12 excluded); the black 
cross marks the accepted value 0.2411 (see text for further details). 



The validity of the methods has been tested with analyti- 
cal calculations on the one dimensional XY model in trans- 
verse field and numerically on a nontrivial spin-1 chain with 
anisotropy. As extreme case, we have shown that the homo- 
geneity condition method provides a satisfactory location of 
the critical point also in the case of Berezinskii-Kosterlitz- 
Thouless transitions. These confirmations motivate us to con- 
sider systems in higher spatial dimensionality. In these cases 
the numerical data are restricted to smaller system-size and the 
need for fast-converging pseudo-critical sequences is a prereq- 
uisite for the precise location of the critical points. 



Appendix A: 

Here we indicate how to compute the mean energy per site 
e^, and the average potential 6^, as require id in Subsec. UlI Al 
Consider for example the energy sum 



L-l 



eL(h,r]) = -— A(fc n ) , k n = 



2n + l 



-7T, 



n=0 



where A(k) is given by Eq. (TT9b . We need to investigate the 
above sum in the quasi-critical region z = L/£ <C 1 so it is 
sufficient to expand ei(/i, rj) in powers of (h — 1) = z/L 



1 L_1 z 

Wjr Y [ A h=l (M + ^ d hMh=l (K, 



n=0 
Z\2 1 



j-) -d 2 h A\ h=1 (k n ) + 0(z 3 ) 



The resulting sums can be computed using the Euler- 
Maclaurin formula (see e.g. lfl4Tl ) 



i=o Ja 



m X 2r 



E 



^(2r)r 



2r 



/< ar - 1 >(&)-/< ap - 1 >(a) 



R m (Al) 



valid for a function / with at least 2m continuous derivatives 
in (a, b). Here B n are the Bernoulli numbers and the remain- 
der R m depends on /( 2m ) on (a, b). Some care must be taken 
when the function / diverges at the border of the Brillouin 
zone, in this case one must keep a and b away from the bor- 
ders. Moreover, sending m to infinity in Eq. dAU . some sums 
must be regularized using a Borel summation technique. The 
final result for the ground state energy is 

fh x n x *Mr-a 7tt 3 (4^-3) r _ 4 
e L (ft, 77) = eoo (1,7?) —L - — .,, L 



+ (h-l) 







boo (M) 



360 4|r7| 



(ft - l) 2 In (L) + 7c + In (8 \rj\ /tt) - 1 
2 n\ V \ 



0(z*) 
(A2) 
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where -fc = 0.577216 ... is the Euler-Mascheroni constant 
and the thermodynamic values are given by 



eoo (M) 



boc{l,v) = — 



2^1 2arctan(yi -if/\r]\ 

TT 



7ryl — r] 2 
axctaa(VWYM; 



iry/1 -rj 2 
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Using similar procedures we obtain the following expression for the evaluation of the finite-size gap in the PRG method, 
for average potential br, (h, rf): can be treated along similar lines. The final result for the gap 



* L -2 7^ 3 (3 - 2 7 2 ) r _ 4 
12 M 2880 |t7| : 



is 



bL (h, i]) = boo (1, rf) - 777777^ 2 , , |3 L 



_ {h _ 1} l„(L) + TC + ln( 8MW -l + a (j2) AL = (h _ 1}+ *(V + j-V, (A4) 
Finally, the sum 192 \r)\ 



L-l 



2nj 



L ^ \ L 

3=0 
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